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We investigate one-dimensional S = 1/2 and 5 = 1 antiferromagnetic quantum spin chains with 
easy-axis anisotropics in a transverse field. We calculate both the uniform magnetization and the 
staggered magnetization, using a variant of the density matrix renormahzation group method. We 
find that the staggered magnetizations exhibit reentrant behavior as functions of the transverse 
field, where the competition between the classical Neel order and the quantum fluctuations plays an 
important role. We also discuss the critical behavior associated with the staggered magnetizations. 
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I. INTRODUCTION 

Recently, magnetization processes of one-dimensional 
(ID) quantum spin systems have been attracting much 
interest from theoretical and experimental view points. 
Particularly, various types of field-induced critical be- 
havior in the magnetization process, such as apScuiarc- 
root behavior associated with the excitationj?apErEl, the 
magnetization plateauQ, the cusp singularity^, etc. have 
been investigated intensively. In most of these theoretical 
analyses of the magnetization process, the external field 
is supposed to be applied in the longitudinal direction. 
For the case of a system having an anisotropy, however, 
magnetic properties of the system may be essentially dif- 
ferent, depending on whether the external field is applied 
in the traverse or longitudinal directions. In this paper, 
we wish to make a quantitative analysis for the magnetic 
properties of antiferromagnetic quantum spin chains with 
easy-axis anisotropics in a transverse field. 

In analyzing an anisotropic model in the transverse 
magnetic field, an essential point is that the "total S z " 
of the system is not a conservative quantity, since the 
symmetry of the model is reduced to be Zi . Then, the 
property of the system is characterized with the staggered 
magnetization. For example, the appearance of the stag- 
gered magnetization is known for the S — 1/2 XXZ spin 
chain in the Ising limit, while the isotropic XXX model 
has no Neel order; The anisotropy favors the classical 
Neel order, whereas the transverse field term induces a 
quantum fluctuation in addition to the usual XY-term. 
What we are interested in here is how the system behaves 
between the Ising limit and the isotropic limit. We dis- 
cuss the physics developed by the competition between 
the anisotropy and the quantum fluctuations, based on 
quantitative calculations of the staggered magnetization. 

For the purpose of computing the magnetization curve, 
we employ the product fiW c function renormahzation 
group(PWFRG) method,&Ll which is a variant of t 
density matrix renormahzation group(DMRG) metho- 



The PWFRG is successfullypapjilied to various types 
of the magnetization process, Lrffl'tHMJ and further it can 
treat a system without the total-,? 2 conservation effi- 
ciently. We calculate the staggered magnetization in the 
easy-axis direction and the uniform magnetization in the 
transverse one, for the infinite-length chains. Here it 
should be noted that the lack of the total-S z conserva- 
tion makes the exact diagonalization study difficult for 
the anisotropic models in the transverse field. 

For the S = 1/2 XXZ chain, we find the reentrant be- 
havior of the staggered magnetization, which is partially 
supported by the exact results at some special param- 
eters. We further discuss the 5 = 1 XXZ spin chain 
and the S = 1 Heisenberg spin chain with the single- 
ion anisotropy, where the existence of the Halda.13 
makes the physics of the systems more fruitful. ElB tiJ 
We verify that the phase transitions occur twice at the 
lower critical field and the higher critical one; The former 
field is associated with the Haldane gap, and the later one 
can be connected with the saturation field of the isotropic 
case adiabatically. We also discuss the critical behavior 
characterized by the staggered magnetizations. 

This paper is organized as follows: In §H, we de- 
scribe the S = 1/2 XXZ spin chain in a transverse field. 
In Wl\ we consider the S = 1 quantum spin chains 



with the XXZ-type anisotropy (^T\- A) and the single- 
ion anisotropy (§[III|-B). The last section is devoted to the 
summary. 



II. THE 5=1/2 CASE 

The S — 1/2 XXZ chain is one of the most fundamental 
models among spin chains of S = 1/2. In this section, we 
consider the XXZ chain in a transverse field. We write 
the Hamiltonian TL as 



qy qv x 



(1) 
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where S" is the a component of the S = 1/2 spin oper- 
ator at i-th site, and T is the strength of the transverse 
field. The exchange coupling constant in easy-plane (xy- 
plane) is denoted by J|| and that of the easy-axis(z-axis) 
direction by J±. We have set g/j,B = 1 (g: ^-factor, /i^: 
Bohr magneton) for simplicity. In this section, we dis- 
cuss the system (Q) in the range < J|| < 1 with fixing 
J± = l. 

The magnetic property_af Ihe Hamiltonian ([!]) has been 
studied for several years .EirEE Particularly, it should be 
remarked that there is a special value of the transverse 
field called a disorder point-j-where the ground state is 
written down analyticallyolij In addition, Mori et al. 
have discussed how the excitation-structure of eq. (||) 
is affected by the transverse field.OEa However, the full 
magnetization process for general parameters of eq.(Q) 
has not been calculated quantitatively, except for the ex- 
actly solved cases Jm = (the Ising model in a transverse 
field) and Jm = 1 (the Heisenberg model). 

In the following, we calculate the uniform magnetiza- 
tion in x-direction M x = (Sf) and the staggered magneti- 
zation in z-direction M st = \{Sf)\, by using the PWFRG 
method. Here we note that it is useful to consider the ir- 
rotation of spins at every two sites(for example, at every 
even sites): 
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since we can calculate the staggered magnetization as 
the uniform magnetization in the transformed system 
Ji! = U'tiU. We have confirmed that the computed re- 
sults for the transformed model agree with the staggered 
magnetization calculated directly for eq. ([!]). 

Before proceeding to the detailed analysis of the XXZ 
chain, we consider the Ising model in a transverse field 
(ITF) ( J|| = 0), for which the analytical solution is known 
as followsEJ: 



M x = 



* dk 



(cos k + 2T) 



and, 



2tt y/l + 4r cos k + 4r 2 ' 
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We can then check the efficiency of the PWFRG method, 
by comparing the obtained results with the above exact 
expressions. In Fig. |l|, we can see that the computed 
results are in good agreement with the exact ones, where 
the PWFRG calculation is performed with the number 
of retained bases m = 14. Thus, we can expect that 
the PWFRG method with a relatively small m yields 
accurate results in the following analysis of the XXZ spin 
chain. 




FIG. 1. The PWFRG results (represented by o) for the 
Ising model in a transverse field: (a) the uniform magnetiza- 
tion M x and (b) the staggered magnetization M s t- The exact 
solutions of M x and M at are denoted by the solid lines. 




FIG. 2. The PWFRG results for the S = 1/2 XXZ spin 
chain with various Ju : (a) the uniform magnetization M x 
and (b) the staggered magnetization M st . The V-M x curves 
correspond to Jm = 0, 0.2, 0.4, 0.6, 0.8, 1 from left to right. 
The r-M st curves correspond to Jy = 0, 0.2, 0.4, 0.6, 0.8, 
from left to right. The symbol o indicates the disorder point. 
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In Fig. g, we display the PWFRG results of the magne- 
tization in the x direction M x for Jy = 0, 0.2, 0.4, 0.6, 0.8, 
1 and the staggered magnetization M st for Jy = 0, 0.2, 
0.4, 0.6, 0.8. In Fig@-(a), the M x grows monotonously 
as the field V is increased, whereas in Fig.^-(b), the 
staggered magnetizations M st for Jy > 0.4 exhibit the 
reentrant behavior. This reentrant behavior can be ex- 
plained qualitatively in terms of the competition between 
the classical Neel order and two kinds of the quantum 
fluctuations. In the low-field region, the system is un- 
der influence of the quantum fluctuation from the XY- 
tcrm( J,, J2 SfSf +1 + SfSf +1 ), while, in the high-field re- 
gion, the effect of the transverse-field becomes dominant. 
Although both of them can decrease the classical Neel 
order, the staggered magnetization is rather enhanced in 
the middle-field region. This enhancement implies that 
the two quantum fluctuations conflict with each other 
rather than reduce the Neel order. This picture is consis- 
tent with the shape changes of the static structure factor 
in the transverse field, as was reported in RefJia. 

Let us now proceed to analysis of details in the mag- 
netization curve. In the zero transverse fieldplipit, the 
staggered magnetization is calculated exactly:E^JE2l 
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where q = 1/ Jy — y 1/ J| — 1. By comparing the PWFRG 

results at T — with the exact values of eq.(||), we have 
confirmed that the obtained data for Jy < 0.6 agree with 
the exact values in 5 digits, within the numbers of the 
retained bases m = 88. 

There is another special value of a transverse field 
called disorder point, wh ere the ground state can be cal- 
culated analytically.ESlia The disorder-point field Td p is 
given by 



dp 



2J||(J|| 



1) 
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and the exact values of the magnetizations are obtained 
as 



(7) 




for Jn > 0.73, by comparing eq.(||) and eq.(||). This 
inequality (0) guarantees the reentrant behavior of the 
staggered magnetization analytically in the range Jy > 
0.73. A precise numerical calculation of Af st (using max- 
imum m = 30) illustrates that the reentrant behavior 
emerges for Jy > 0.3984. 

We next discuss the critical behavior of eq. ([[]) . The 
critical field T s is defined by the field at which the stag- 
gered magnetization vanishes. Since the symmetry of the 
system is Z2, the universality class is expected to be the 
two-dimensional(2D) classical Ising one, where the stag- 
gered magnetization behaves as 



Ms, 



(r s -r)' 3 for r<r s 



(10) 



with (3 — 1/8. For the case of Jn = 0, we can see (3 = 
1/8 from the exact solution (f|) easily. For < Jy < 
1, we have checked the eq.([l0]) with the PWFRG result 
of the staggered magnetization. For example , we show 
the plot of r-(M st ) 8 for Jy = 0.8 in Fig]|, where we 
can see the good linearity of the plotted data. We have 
performed the fitting of M st — A(T S — T) 13 with fitting 
parameters A, (3 and T s , and then obtain the exponent 
j3 = 0.124,0.124 and 0.125, for Jy = 0.4,0.6 and 0.8, 
respectively. These values are consistent with (3 = 1/8. 
Thus, we have verified that the universality class of the 
XXZ model at the critical transverse field is the 2D Ising 
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FIG. 3. The check of linearity of the r-(M st ) 8 curve in the 
case of Ju = 0.8 near the critical field. The good linearity(see 
eq. (jlo|)) can be seen. 
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III. S = 1 CASES 



Af-t = r< 



l l-Jn 



1 + Jn 



(8) 



In Figs §-(a) and (b) we have plotted the disorder point 
with the open circles, which are consistent with the 
curves calculated by the PWFRG. 

Here we remark that we can obtain an inequality: 



In this section, we consider the 5 = 1 XXZ spin chain 
and the 5=1 Heisenberg spin chain with a single-ion 
anisotropy. In contrast to the 5 = 1/2 XXZ model, the 
5 == 1. ch ains with small anisotropics have the Haldane 
gapEraOEj. Then it is an interesting problem to clarify 
how the Haldane gap influences the competition between 
the classical Neel order and the transverse field. 
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A. the S = 1 XXZ model 

In this subsection, we consider the S = 1 XXZ spin 
chain, whose Hamiltonian is written as 

w = EWW+i + + J±s?s? +1 ] - r^sf, 

(ii) 

where iS" is the a component of the 5* = 1 spin operator 
at i-th site. In this section, we fix Ju — 1 and vary J± in 
J_l > 1, for the later convenience. 




FIG. 4. The PWFRG results for the S = 1 XXZ spin chain: 
(a) the uniform magnetization M x and (b) the staggered mag- 
netization M st . Each curve corresponds to J± = 1.5, 1.4, 1.3, 
1.2, 1.15, 1.1, and 1.05 from top to bottom. The open circle 
represents the disorder point for each J±. 

In Fig.^J, we show the staggered magnetization M st and 
the uniform magnetization M x . The maximum number 
of the retained bases in the DMRG calculation is m = 70. 



In the figures, we . 
S = 1 XXZ model,! 



sa plot the disorder points for the 
EL3 which is given by 



M x = 



1 + J A 



and M a t = 



J 



1 



1' 
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with the disorder-point field Tdp = 2^/2(1 + J±). We 
can see that they are in good agreement with the DMRGj 
results. When the anisotropy is not so big(Ji < 1.18),c3 
the Haldane state is maintained up to the critical value of 
the transverse field r c i , where the phase transition occurs 
accompanying the spontaneous staggered order(Fig.4- 
(b)). As r is increased above r c i, the staggered magneti- 
zation grows rapidly. However, as T is further increased, 



M s t turns to decrease and finally vanishes at the higher 
critical field T C 2- In Fig. f§-(a), we show the uniform mag- 
netization Mr, where the anisotropy effect seems not to 
be significant. However, a notable point for this uniform 
magnetization is that M x has a small but finite value be- 
low the critical field r c i , in contrast to the isotropic case 
where M x is exactly zero below r c i . In a real experi- 
ment, it may be important to distinguish such a behavior 
of M x from the finite temperature effect. 

As J± is increased, the lower-critical field r c i be- 
comes small and thus the region of the Haldane phase 
shrinks (see FigJ^-(a)). When J± is increased above the 
critical value of the anisotropy J± c ~ 1.18, the sponta- 
neous staggered magnetization appears at the zero trans- 
verse field£3 and M st ~ ( J± — J± c ) 1 ^ 8 is observedo 

For J± > Jie the staggered magnetization exhibits the 
week reentrant behavior, which is similar to the S = 1/2 
case. In the J± — > oo limit, the system becomes to the 
5=1 version of the Ising model in a transverse field, 
without the reentrance of M st . 
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FIG. 5. Critical behavior of the staggered magnetizations 
for the 5" = 1 XXZ model: (a)S = 1 XXZ chain with J± = 1.2 
near the higher critical field V C 2 and (b) 5* = 1 XXZ chain with 
J± = 1.1 near the lower critical field T&. 

Let us next discuss the critical behavior at T = T c i and 
r c 2, which are also expected to belong to the 2D Ising 
universality class. In Figs.||-(a) and (b), we show the plot 
of (M st ) 8 near the critical fields r c i and T C 2, respectively. 
In these figures, we can recognize the good linearity. Near 
the lower critical field r c i , we have estimated the critical 
exponent f3 = 0.126 for J± = 1.1. Also near the higher 



4 



critical field T C 2, we have (3 = 0.126 and 0.127 for J± = 
1.2 and 1.4, respectively. These values of (3 are consistent 
with the 2D Ising universality class. The critical behavior 
of M x should correspond to that of the internal energy 
of the 2D Ising model, as well. However it is difficult to 
estimate the exponent a with a direct calculation, since 
the expected value a = 0(log) is too weak to be detected 
numerically. 



B. S = 1 Heisenberg model with a single-ion 
anisotropy 

For the case of the 5=1 quantum spin chains, wc 
can consider the single-ion-type crystal field as an origin 
of the anisotsppy, which is actually observed in the real 
compounds.E-3 The Hamiltonian is given by 



h = £ 3 • 3+i + D £( 5 ?) 2 - r £ s i 



(13) 



If D < 0, the £)-term( J DX; i (5 , f) 2 ) stabilizes "±1" spins 
equivalently and hence the D-term is expected to make 
the similar effect to the previous XXZ-type anisotropy. 
In the following, we thus consider D < region, where 
the staggered order in the z-direction is induced. The 
phase boundary-, of this model is investigated by Sakai 
and Takahashili3, using the exact diagonalization up to 
14 sites. However the magnetization curve of this model 
has not been calculated yet. 




FIG. 6. The PWFRG results for the S = 1 Heisenberg 
model with the single-ion anisotropy: (a) the uniform mag- 
netization M x and (b) the staggered magnetization M s t- In 
Fig. (a), the curves correspond to D — —0.4, —0.3, —0.2, 
—0.1, and —0.05 from left to right. In Fig. (b), the curves 
correspond to D = -0.45, -0.4, -0.35, -0.3, -0.25, -0.2, 
—0.15, —0.1, and —0.05 from top to bottom. 

In Fig.[| we show the uniform magnetization M x and 
the staggered magnetization M st . We can see that they 
exhibit the similar behavior to that of the S = 1 XXZ 
models. For —0.3 < D < 0, the system has the staggered 
order between the two critical fields r c i and T C 2- In T < 
r cl , the system is in the Haldane phase, but the uniform 
magnetization M x has a finite value due to the lack of 
the total-S z conservation law. As D is decreased, r cl 
becomes small, and finally vanishes at the critical value 
of the anisotropy D c ~ —0.3, which-ds consistent with 
D c = —0.29 calculated by Sakai et al.l!3 As D is decreased 
further, the staggered magnetization appears at the zero 
transverse field, and then the reentrant behavior of the 
staggered magnetization can be seen. 

Next, we plot (M st ) 8 near the critical fields in Fig. [?], 
where the good linearity can be seen. For example, the 
estimated critical exponent at the lower-critical field is 
f3 = 0.127 for D = —0.15, and the ones at the higher 
critical fields are j3 = 0.127 and 0.126 for D = -0.15 and 
—0.45 respectively. These values are also consistent with 
the universality class of the 2D Ising model. 
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FIG. 7. Critical behavior of the staggered magnetizations 
for S = 1 Heisenberg model with the single-ion anisotropy of 
D = —0.15: (a)near the higher critical field F C 2 and (b) near 
the lower critical field r c i. 



IV. SUMMARY 

In this paper, we have considered the one-dimensional 
quantum spin chains with the easy-axis anisotropy in the 
transverse magnetic field. We have calculated the stag- 
gered magnetization M st in the z direction and the uni- 
form magnetization M x along the x direction precisely, 
using a variant of the DMRG(thc PWFRG method). For 
the 5 = 1/2 XXZ chain, we have shown that the M st ex- 
hibits the reentrant behavior when Jm > 0.3984, where 
the competition between the classical Neel order and the 
quantum fluctuations from the transverse field and the 
YY-tcrm plays an important role. We have also dis- 
cussed the critical behavior associated with the staggered 
magnetization and verified that the universality class of 
the model is the two-dimensional(2D) Ising one. 

For the 5 — 1 XXZ model and the 5=1 Heisenberg 
model with a single- ion anisotropy, we have shown that, 
if the anisotropies are not big(Ji < 1.18 or D > —0.30), 
the systems have the staggered order between the lower 
critical field associated with the Haldane gap(r c i) and 
the higher one connected to the saturation field of the 
isotropic case(r c 2). When the anisotropies are big, the 
staggered order appears at the zero transverse field with 
the reentrant behavior, which is similar to that of the 5 = 
1/2 XXZ model. We have also investigated the critical 
behavior at r c i and T C 2, by estimating the exponent (5. 
As a result, we have verified that they are consistent with 
the class of the 2D classical Ising model. 

In the actual one-dimensional compounds, various 
types of the anisotropy can be considered. When ana- 
lyzing the magnetization processes of such matters, we 
believe that the present results are useful. 
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